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Abstract 

We study the spectral measure of large Euclidean random matrices. The entries of these 
matrices are determined by the relative position of n random points in a compact set fin 
ofR^. Under various assumptions we establish the almost sure convergence of the limiting 
spectral measure as the number of points goes to infinity. The moments of the limiting 
distribution are computed, and we prove that the limit of this limiting distribution as the 
density of points goes to infinity has a nice expression. We apply our results to the adjacency 
matrix of the geometric graph. 

Keywords: random matrix, spectral measure, random geometric graphs, spatial point 
process, Euclidean distance matrix. 

1 Introduction 

The main research effort in the theory of random matrices concerns matrices where the coeffi- 
cients are independent random variables (see Bai [Ij for a survey). Few authors have studied 
the limiting spectral measures of other types of large matrices, in particular, Markov, Hankel 
and Toeplitz matrices have been studied by Bryc, Dembo and Jiang [3] and Toeplitz matrices 
by Hammond and Miller |10]- In this paper, we consider another class of random matrices, the 
Euclidean random matrices (ERM) which have been introduced by Mezard, Parisi and Zee |13j . 
An ERM is an n X n matrix. A, whose entries is a function of the positions of n random points in 
a compact set Q, of M'^. In this paper, O will be an hypercube, the n points Xn = {^i, • • • , ^n}, 
n uniformly distributed points in J7 and 

A = {F{Xi — Xj))i<i<j<n, (1) 

where F is a measurable mapping from to C. We will pay attention to the spectral properties 
of A. In this paper, we will compute some limits of the spectral measure as the number of points 
n goes to infinity. We will show how the eigenvalues of A are related to the Fourier transform 
of the mapping F. 

Examples of interests in branches of physics are explained in [13] and Offer and Simons [H] . 
A particularly appealing case is F{x) = \\x\\, the Euclidean norm. This subclass of ERM is 
called the random Euclidean Distance Matrices and some of their spectral properties are derived 
in Vershik [17] , Bogomolny, Bohigas and Schmidt [2] . In pPT] the author generalizes the problem 
considered here and consider an integral operator defined on a metric space. Koltchinskii and 
Gine [Tl] have analyzed the convergence of the spectra of the matrix (n^^/i(Xj, Xj))i<jj<„ to 
the spectra of the compact integral operator with symmetric kernel h. 
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Another field of application is graph theory. Indeed, if F{x) = 1(0 < ||x|| < r), then A 
is the adjacency matrix of the proximity (or geometric) graph (refer to Penrose [15j). More 
generally if F{X) = F{—X) G {0, 1} then A is the adjacency matrix of a random graph. The 
spectral properties of the adjacency matrix or related matrices are of prime interest in graph 
theory. For example the probability of hitting times of random walks on graphs is governed by 
the spectrum of the transition matrix (for a survey on this subject, see e.g. Section 3 in Lovasz 
[12|). Or, in network epidemics, the time evolution of the infected population is also closely 
related to the spectral radius and the spectral gap of the adjacency matrix, see Draief, Ganesh 
and Massoulie [8]. For Erdos-Renyi random graphs, some properties of the spectrum can been 
computed thanks to the seminal work Wigner of [18] and Fiiredi and Komlos [9]. For power 
law graphs and related graphs, see Chung, Lu and Vu [1], [5]. 

Various generalizations of ([T]) would be worth to consider. Some extra randomness in the 
model could be added, and the entry of the matrix i,j could be equal to Fij{Xi — Xj), where 
{Fij)i<i,j<n are i.i.d. mappings independent of the point set X^- Falls into this framework 
the adjacency matrix of a random graph where there is an edge between two points with a 
probability which is deterministic function of their distance, such as the small world graphs 
(see for example Ganesh and Draief [7]). 

Another generalization is the original model of Mezard, Parisi and Zee [13] where the entry 
i,j is equal to 

F{Xi - X,) - u5,j Pi^i - ^k), 

k 

where 6ij is the Kronecker symbol and u G M. The case u = 1 is of particular interest, the 
matrix is then a Markov matrix. 

In order to obtain the adjacency matrix of more sophisticated geometric graphs, such as the 
Delaunay triangulation, it would be necessary to consider an entry i,j which depends on the 
whole point set Xn and not only on Xi — Xj. 

We will consider two models in this note. In the first model, Q. = [—1/2, 1/2]'^ and F is 
1-periodic function: if x,y E M'^ and x — y € then F{x) = F{y). Equivalently, the point 
set Xn = {Xi, • • • , Xn} could be on the unit torus T'' = W^\Z'^. We choose a periodic function 
in order to avoid all boundary effects with the hypercube 0. The matrix A is defined by 
where -F is a measurable function from M'^ to C. The discrete Fourier transform of F is defined 
for all A; E Z*^ by F{k) = F{x)e~'^™^'^dx. Throughout the paper, we assume that a.e. and 
at 0, the Fourier series of F is equal to F: 

F{x) = F{k)e^'''^-''. 

keZ'' 

A sufficient condition is YlkeZ'^ 1-^(^)1 < and F continuous at 0. This Fourier transform 
plays an important role in the spectrum of A. As an example, consider U = {Ui)i<_i<n a vector 
m and assume F hermitian {F{—x) = F{x)), then a.s. 

" 2 

U*AU = - Xj)UiUj = YY1 F{k)e^'"''-^^'~^'^U^Uj = Y ^i^) \ Y e^'"'''^'^^ | ■ 

i,i i,j k k 1=1 

Therefore A is positive if and only if for all k E Z'^, F{k) > 0. 

We will compute explicitly the spectral measure of the matrix An = A/n as n tends to oo, 

n 
i=l 
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where {Aj(n)}i<j<„ is the set of eigenvalues of A. Notice that {Aj(n)/n}i<j<„ is the set of 
eigenvalues of An- We define the measure: 

Since limn^jn^o^^ F'ik) = 0, ^ is a counting measure with an accumulation point at 0. 
Theorem 1 For all Borel sets K with fi{dK) = and ^ K, a.s. 

limf,n{K) = fi{K). (2) 

n 

The convergence of the spectral measure follows also from Theorem 3.1 in [TT]. As an 
immediate corollary, we obtain the convergence of the spectral radius of An, almost surely, 
limji^oo maxi<j<„ ^^^^ = maxj^^^d \F{k)\. For example if F{x) = l(maxi<j<rf < r) then 
F{k) = r'^Y[i=i sinc{27Tkjr), where sinc(x) = sin(x)/x and k = {ki, ■ ■ ■ , kd) G Z*^- The spectral 
radius of An converges a.s. to r'^ and the second largest eigenvalue to r°'sinc(27rr) if r is small 
enough, thus the spectral gap is equivalent to r'^'^^(27r)^/3! as r goes to 0. 

Our second model is more challenging. Now, Xn = {Xi, ■ ■ ■ , Xn} is the set of n independent 
points uniformly distributed on the hypercube S'^Q = [—5~^/2, /2]'^ where 5„ goes to 0. In 
this second model, we scale jointly the number of points and the space. We assume that for 
some 7 > 0, 

limJ^n = 7. (3) 

7 is the asymptotic density of the point set Let / be a measurable function from to C 
with support included in 0, the matrix A is defined by ([T|) (with F replaced by /). 

Considering the change of variable x 5x, the matrix A is equal to the matrix Bn defined 

by 

Bn = ifSniXi — Xj))i<i<j<n, 

where fs'-x^^ f{^/^) and the point set Xn = {-^i, • • • iXn} is a set of n independent points 
uniformly distributed on Q.. The spectrum of Bn is denoted by (A']^(n), • • • , A^(n)), we define 
the empirical measure of its eigenvalues: 

1 " 

= - / ,^\'-(n)^ 
i=l 

We will prove the following: 

Theorem 2 For all ^ > 0, there exists a measure such that for the topology of the weak 
convergence, a.s.: 

lim Un = l^-y- 

Moreover 71-^1/^ is continuous (for the topology of the weak convergence) . 

The exact computation of is a difficult problem, we will compute the value i^^(Pm); 
where Pm is the polynomial t ^ (Equation ()22p l. However, the behavior of as 7 goes 
to infinity is simpler. Indeed, we define the Fourier transform of / by, for all ^ E R'^, f{C) = 

e~'^™^'^ f{x)dx. Since / has a bounded support, / is infinitely differentiable. We assume 
that the following inversion formula holds 

a.e. and at 0, f{x) = [ /(Oe^^^^'^''^^ (4) 
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Note that if / is hermitian (/(-x) = f{x)) then /(^ G K and for e > 0, J^^ > e)dC 

is finite. Hence by the change of variable formula, there exists a function ■0 such that for all 
continuous functions h with ^ supp{h): 

[ h{t)i){t)dt= [ h{f{x))dx. 

ip is the level sets function of /, if £ denotes the d-dimensional Lebesgue measure, for all t > 0, 
ip{t) = lim^^Q £{{x : \f{x) — t\ < e})/e. If c? = 1 and / is a diffeomorphism from M to -ftT then 
Tp has support on K and is equal to ip{t) = 

Theorem 3 If f is hermitian and holds true, then as 7 goes to infinity, for all analytic 
functions h{t) = X^meN ^m*™ ^^^^ ^0 = and Ylrnen l^ml^™" finite for all t: 

[ h{t)u^{dt) ~ f h{t)-f-'^ij{-)dt = [ -f-^h{jt)ij{t)dt. 
Jr Jr 7 Jr 

This result states that the measure v^{dt) is in a weak sense equivalent (not in the measure 
theory sense) to the measure 'y~'^ip{t/j)dt in the high density asymptotic. The measure ip{t)dt 
is the continuous analog of the counting measure n in Theorem [TJ As an example, if d = 1 and 
f{x) = 1(0 < |x| < r) then /(^) = rsinc(27r^r) and 'ijj{t) is plotted in Figured) 




' O 2 4 6 8 10 -0.2 O 0.2 0.4 0.6 0.8 



Figure 1: Left: /(^) for f(x) = 1(0 < |x| < r). Right: the level set function ip. 



Remark. Let 7^ = n6'^, if 7^ tends to infinity and 5„ goes to 0, with the material of this note, 
we may also prove the convergence of Yl^=i ^x'{n)/jn *° measure '4){t)dt on all continuous 
function h with compact support and ^ supp{h). 

The spectral radius of the matrix Bn is not computed explicitly in this paper. However, the 
following upper bound is available: 

Proposition 4 If d > 2 and Po{^) denotes a random variable with Poisson distribution of 
intensity 7, then with a probability tending to 1 as n goes to infinity, 

max |A-(n)| < j(?^) sup |/(x)|, 
where j{n) is solution of: nP{Po{'~f) > j{n) + 1) < 1 < nF{Po{'~f) > j{n)). 
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For n large enough, using the inequahty F{Po{'y) > k) < exp(— | ln(^)), for k > e^7, we deduce 
that j(n) < 3 Inn/ In Inn for n large enough. 

The remainder of this paper is organized as follows. In Section [2l we prove Theorem [H 
In Section [3l we prove Theorems [21 [3] and Proposition [H Finally, in Section HI we state some 
simple results on the eigenvectors of A and on the correlation of the eigenvalues. 

By convention C will denote a constant which does not depend on n. Its exact value may 
change throughout the paper. Also we define: ||-F||oo = sup^.g]gd and B{x,r) will denote 

the open ball of radius r and center x on the torus T*^. 



2 Proof of Theorem [T] 

The proof of Theorem [T] relies on the classical Wigner's method [TS] to compute the empirical 
mean distribution measure of eigenvalues. We will compute for all m G N: 



1 1 " 

lEtr^n = -^^trA- = — ^ = f^niPm). 



n n 

i=l 



We will then use a Talagrand's concentration inequality to prove that trA"^ is not far from its 
mean and conclude. About the rate of convergence of ^u^ to fi, we will state (in the forthcoming 
LemmaE]) that, if Pmit) =t'^,m> 1, 

m—i rn(m 1) 

limn ( EfXniPm) - /U(P„^) )= qiliPg)fliPm-q) - V (^m)- (5) 

q=l 

We begin with a technical lemma. 

Lemma 5 For < p < m, let Sm,p be the set of surjective mappings from {!,••• ,m} to 
{1, • • • ,p}. We have: 

m ^ \ „ m 

^trA"^ = ^r) E / ll^(^m-^Hj+i))dxi---dx„ (6) 

with (j){m + 1) = 0(1) and with the convention that (^ = for p > n. 
Proof. By definition: 

m 

ii,---im j=l 

with im+i = ii and the sum is over all n-tuples of integers i = (ii, • • • , im) in {1, n}*". Let p{i) 
be the set of distinct indices in i. We can define a surjective mapping (pi in such that 

ij = i(j,-^[j)- Taking the expectation in Equation we get 

„ m 

EtrA™ = Y Ii ^(^Mj) ~ ^</.i(j+i))'^2;i • • • rfxp(i), 

i=(ii,---jm) ^ i=l 

We then reorder the terms. We consider the equivalence relation in ^ 0' i-^^ there 

exists a permutation a of {1, • • • ,p} such that a o cp = cp'. The value of J^^ TYJLi ^(^<t>(j) ~ 
x^(^j^i))dxi ■ ■ ■ dxp is constant on each equivalence class. Let (j) S the numbers of indices i 

such that (pi ^ (p is equal to n!/(n — p)! (if n > p and otherwise). Since there are p\ surjective 
mappings in the class of equivalence of (p, we deduce Equation □ 
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Lemma 6 For each m, 



EyUn(Pm) = + ' ( Yl Ql^iPq)KPm-q) " ^> (Pm) ) +o(-) 

Tt \ ^ /in 



Proof. We apply Lemma [S] and identify the coefficients in and n™ ^ in Equation (0) . We 
first consider the term in n"^, such a term comes from p = m: 



/ I I F{xj — Xj+i)dxi • • • dxr^ 
Jn"^ 111 



(n - m)! Jn-^fj{ 
By induction, we easily obtain that 

/ TT ^(^i - Xj+i)dxi ■■■dxm= / F*'^(0)dxi = F*'"(0), 
in™ in 

where * denotes the convolution operator: F*G{y) = /q F{y — x)G{x)dy and F*"^ is F^^i^ • • •*F 

(m times). We recall the two properties: F*G{k) = F{k)G{k) and F(0) = J2keZ'' ^i^)^ ™ 
order to get: 

~ m 

We thus deduce that: 

limE/i„(Pm) = n{Pm)- 

n 

It remains to identify the terms in n™""^ in Equation ([6]). This term comes from two contribu- 
tions p = m and p = m — 1. Since n\/(n — m)\ = n™" — n*""^ SI^o^ ^ "I" o{n'^~^), the term in 
jfn-i in p = m is equal to: 

- ^^^MP™). (8) 

The leading term for p = m — 1 is 



n! 

(n — m + 



Now if G Em,m-i! </> "^{i) is not reduced to a single point for a unique index i^. Since 
the value of /^^m-i HJli Pi^<l>U) ~ ^<t>{j+i))dxi ■ ■ ■ dxm-i is invariant under permutations of the 
indices, without loss of generality, we may assume that = 1, = + 1} with 

(7 G {1, • • • , m — 1} and </>(j) = j if j < (7 and = j -)- 1 if j > g + 1. For such 0, integrating 
over X2, • • • • • • ,Xm-i, 



/ IT P{x<l>{j) - X^fj+i))dxi ■ ■ ■ dXm- 



i7*('?)(0)F*('"-«)(0)dxi 



Finally, for each q, there are {m—q) x (m— 1)! surjective mappings such that, up to a permutation 
of the indices, (t)~^{l) = {1, g + 1} and = j if j < g and = j + \ li j > q + 1. Indeed for 
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such 0, there are (m — q) possible pairs (ii, ii + g), 1 < ii < m — q such that (pih) = (p{ii + 
Therefore Equation ([9]) can be written as: 



n!(m — 1)! 



m—l 



m—1 



(n — m + l)!(m — 1) 



— ^{m- q)fi{Pq)n{Pr 



m~q j 



n 



q=l 



''Y.qf,{PMPm-q) + o{n"^-'). (10) 

9=1 



Adding this last term with the term ([8]), we get the stated formula. 

We may now prove Theorem [TJ 
Proof of Theorem Ul 

We fix n and for each m > 1, we define the functional: 



□ 



Qm i'^n) 



1 



n 



m— 1 



tr^'" = nflniPn 



If X, y G fi", let (i(x, y) = Yli=i ^{^i Hi) denote the Hamming distance. The functional Qm 
is Lipschitz for the Hamming distance d. Indeed, define x' = (x^)i<j<„ by x^- = Xj for j ^ I 
and x\ ^ Xi, we have: 



(3m(x) - Qm(x^) 



n 



ii,---,imj=i i=l 



< 2m||/||-, 



indeed | Y[T=i ^i^i-, ~ ^ii+i) ~ Y[T=i ^i^i ~ ^i -^i )l most 2||/||™ and it is non zero only if 
there exists a index ij such that ij = I- It follows easily that Qm is 2m||/||™ -Lipschitz for the 
Hamming distance d. 

Let Mm denote the median of Qm- We may apply a Talagrand's Concentration Inequality 
(see for example Proposition 2.1 of Talagrand |16j). 



Fi\Qm-Mm\ >t) <4exp(- 



4m^ 



2m 



n 



integrating over all t we deduce: 

\nEfi{Pm) - Mm\ < n^Pm) " M„| < CmV^, 

for some constant Cm and it follows, that for all s > Cm/V^: 

{S-Cm/V^f. 



\Hn{Pm) - E/i(Pm)| > s) < 4exp(-n- 



4m2||f||2™ 



Using the Borel Cantelli Lemma and Lemma[6l a.s. lim„/i„(Pm) = fi{Pm)- 



□ 



3 Limit Spectral Measure of Scaled ERM 
3.1 Proof of Theorem [2] 

The study of the first model was simplified by the absence of boundary effects with 17. So in 
order to prove Theorem [21 we will first discard them in the second model. We define Fg as the 
1-periodic extension of fs'- for all x G M'*, there exists a unique couple (y, u) such that x = y + u, 
with M G Z'^ and y G 0, and we set Fs{x) = fs{y). 
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We now introduce a matrix and its spectral empirical measure: 

1 

Bn = {F5„ {Xi - Xj))i<i<j<n and = - ^^J 



i=l 

where (Ai(n),-- - ,A„(n)) is the spectrum of Bn. The next lemma states that the limiting 
spectral measures of P„ and Vn are equal. 

Lemma 7 For the topology of the weak convergence of (signed) measures, a.s. z^„ — P„ converges 
as n goes to infinity to the null measure. 

Proof. It is sufficient to prove that for all m > 1, a.s. lim„z/„(Pm) — i'n{Pm) = 0. To this end, 
we notice that if x,y E il, fs{x — y) = Fs{x — y) unless x £ — S)Q and y £ B{x,6). We 
write: 



m m 

l^n{Pm)-i>n{Pm)\ < " E I 11 Z^!! " ^^+1 ) " 11 " ) I 

j=l 3=1 
^ m 

- 2\\f\\Zl{X,,£il\{l-6n)^)\{l{X,^£B{X,,,m5n)) 



< 

n 

< -\\f\\zNnmi-m5n)nr, 



where Nn is the counting measure Nn{-) = G {!)■■■ )^} '■ Xi G ■}. Note that P(Xi G 
fl\{l—m6n)^) < C6n. By the strong law of large numbers, it follows easily that that Nn{0,\{l — 
mSn)^)/n converges almost surely to 0. □ 

By Lemma[71 we may focus on Bn and z>„. In order to keep the notations as light as possible 
we drop the " ~ " in Bn and Dn- 

We first prove that, 

f„ converges in probability to a measure ly-y for the weak convergence. (11) 
By Lemma El if m > 1 , 

1 ^rl\ f 

El^n(-Pm) = -E ( J E / IlPsAx^(j)-X^^j+i))dxi---dXp. (12) 

We begin with an elementary lemma. 
Lemma 8 If (p £ S^.p, p > 1 the value of 



r. m 

/ IT Pi^<p(j) - x^(j+i))dx2 ■■■dxj: 
Jnp~^ ,=1 



does not depend on xi. 



Proof. We consider the change of variable, for j > 1, Xj — Xj X\. The Jacobian of this 
change of variable is 1. If we set x'^ = 0, we obtain /^p-i HJli Pi^HJ) ~ ^<t>{j+i))d'X2 • • • dxp = 
Inp-^ Hj=i ^^^'<t>(j) ~ ^0(j+i))'^^2 ■ ■ ■ ^^p- '-' 
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Assume m > 2, by Lemma [HI we have: 

1 f ri\ f 

EMPm) = FsM'^ + 'Y^ij Yl / llFsAH{j)-H{j+i))dx2---dxp 

1 / \ I, m 

"p=2 V^'^^es™,^ j=i 

where A((/.) = Jljli - x^q+i)) - njLi fSni^HJ) " 2;^(j+i))dx2 • • • (ixp. Since the 

support of fs is included in 6n, if fsix^i^j) - x^(j+i)) 7^ ^<5(^0(j) " ^^^n G 

(5)0. Moreover notice that if Hjli -^(^<A{i)~^0(i+i)) ^ then X2, • • • Xp S B{xi, {m—l)5). 
By Lemma [HJ from now on, we can assume without loss of generality: 

xi = 0, 

and then A((^) =0 for 5 < l/(2m). 

Considering the change of variable yi = Xi/5n in the integrands of Equation (|13p . we obtain, 
for 6 < l/(2m), with yi = 0, 

Finally, since (p) ~ n^/pl as n goes to infinity, we deduce that, for m > 2, 
hm Ei/„(P„) = f{Or + E ^ E / n ^(2^</'(i) - yHJ+i))dy2 ■ ■ ■ dyp. (15) 

(For m < 1, we have z^n(-Po) = 1 and Vn{Pi) = /(O))- 

We check easily that the right hand side of Equation (|15p is bounded by (Cm)™ for some 
constant C not depending on m. Therefore, by Carleman's Condition, there is exists a unique 
measure such that lim„^oo lE^'„(Pm) = i^-y{Pm)- In particular, the sequence (t'n)neN is tight 
and we have proved (jlip . 

The continuity of 7 i— > z^^ follows from the comtinuity of 7 ^ ^^{Pm)- Indeed, let (7n)neN 
be a sequence converging to 7 < 00. Since sup„i^^„(P2) < 00, the sequence (i^7„)n6N is tight. 
Hence for all e > 0, there exists a compact set K such that for all n v^^[K'^) < e. Now, let h 
be a continuous function with compact support, we need to prove that lim„^oo z^7„(^) = i^-y{h). 
Fix e, there exists a polynomial P such that sup^g^^ \ h{x) — P{x)\ < e, we deduce that |z/7„(/i) — 
^^7(^)1 < WyAf^)-i^^AP)\ + W^AP)-^^iP)\ + \MP)-Mh)\ < 2e{l + \\h\U + \u^,SP)-^yiP)\- 
Letting n tends to infinity, since e is arbitrary small and 7 1-^ T^'y{P) is continuous, we obtain: 
lim^^oo z^7„(/i) = v^{h). 

It remains to prove the almost sure convergence of We will prove that for each m > 1, 
there exists a constant C and 

E ( trB;^^ - EtrS™ )^< Cn^ . (16) 

This last equation implies E ( Vn{Pm) — Ei^n(-Pm) C /in? and by Borel Cantelli Lemma, we 
deduce that Vn{Pm) converges almost surely toward u^{Pm)- 
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It remains to prove Inequality (|16p . A circuit in {1, • • • ,n} of length m is a mapping 
TT : Z ^ {1, • • • n} such that for all integer r, 7r(m + r) = 7r(r). Following Bryc, Dembo and 
Jiang [3J, we introduce the new notation: 



Ftt — -^<5„ (-''^7r(i) - ^7r(i+l))- 
1=1 

We then write: 

4 

E(trB--Etri?-)'=E( J^F^-EF^)'= E [Hf^^ -EF^,], (17) 



where the sums are over all circuits in {1, • • • , n} of length m. 

Notice that E [ nf=i -^ti ~ ^^ni ] = if there exists a circuit vTfc, 1 < fc < 4 such that the 
image of vTfc has an empty intersection with the union of the images oi ni, I ^ k. Indeed, due to 
the independence of the variables (Xj)i<j<„, Ft^j. — EF^rj. is then independent of Yli^f, Fn-, —EFt^^. 

Two circuits vri and tt2 in {1, • • • ,n} of length mi and m2 with a non empty intersection 
of their images may be concatenated into a circuit in {1, • • • , n} of length mi + m2 as follows. 
Assume that 7ri(io) = 7r2(jo)) we define the circuit 7ri^2 of length mi+m2 by for i £ {!,■ ■ ■ ,mi + 
m2} 

7ri(i) if I < i < io 

7ri.7r2(i) = ■( vr2(z-io+jo) if io + 1 < « < «o + "^2 

7ri(i — m — io) if io + m+ l<i< mi + m2 



We have: 

Using the same reasoning as for Equation (fT^ . we get 



F F = F 

TTi 7r2 ^ 7ri-7r2* 



» 2m 

EF^^F^2 = 5^^^"^) / TT /(yTTi.TraO) - y7ri.7r2(j + l))%2 " " " > 

where is g = q{iTi,TT2) is the cardinal of the union of the images of vri and 112 and {yi^ , • • • , y^,) 
is the image of vri " 7r2 and yi^ = 0. 

If A^(7ri,7r2) is the cardinal of the intersection of the images of vri and 712, if A^(7ri,7r2) > 1, 
we obtain 

Otherwise, A^(7ri,7r2) = 0, if (/(vTj) is the cardinal of of the image of tTj, 

E\F^^F^^\ = E|F^JE|F^2| < Cn-3(^i)-9(^2)+2 ^ (<.^-g(^i,^2)+2^ ^^^g) 

Similarly assume that A^(7ri, 7r2) > 1, A^(7ri.7r2, tt^) > 1, if (/(vri, 7r2, tt^) is the cardinal of the 
union if the images of 711,-^2, its then F^^^Ft^^F^^.^ = F(^^^ .^^y^^ and we deduce similarly 

E\F^,F^.,F^.^\ < C7n~5(^i'^2,^3)+i_ 

Finally assume that A^(7ri,7r2) > 1, A^(7ri.7r2, tts) > 1, A^(7ri.7r2.7r3, 714) > 1, if (/(vri, 7r2, vra) is the 
cardinal of the union if the images of tti , 7r2 , vrs , 7r4 , we obtain: 

ElF^.F^jF^gF^J < C7n-'?(^i'^2,^3,^4)+i_ (20) 
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By p!7|) . it remains to decompose: 

4 

4! Yl ^[ll^.,-^F^,] 

(7ri,--,7r4)GSUS' 1=1 

where S is the set of quadruples of circuits such that A^(7ri,7r2) > 1, A^(7ri.7r2, vra) > 1, 
A^(7ri.7r2.7r3, 7r4) > 1 and 5' is the set of quadruples of circuits such that A^(7ri,7r2) > 1 and 
A^(7r3,7r4) > 1 and otherwise for i < j, N{ni,nj) = 0. 

The decomposition of the E [ nf=i ^-^i ~ ^^ni ] gives rise to four types of terms: 

2- I](7ri,--,7r4)e5US'^nf=l^vrp 

4- E(,ri,- ,7r4)esus' ^^^h ^^''h ^-^h ^-^lA ' 

where (Zi, ^2, ^3, ^4) is a permutation of (1, 2, 3, 4). We will apply successively the same method 
to bound these terms. 

We begin with the terms of type 1, we have: nf=i ^-^t, < Cn~^'=i'^^'^'^"'"^. Since (tti ,•• • ,114) G 
S U S", Q'(vri,7r2,7r3,7r4) < Ya=i li'^i) - 2, hence: 

4 

JJEF^, < C'n-9('riA2,7r3,7r4)+2^ 
1=1 

There are at most Cn'^ quadruples of circuits such that q'('/ri, 7r2, vrs, 7r4) = q, therefore the terms 
of type 1 may be bounded as by 

4 

(tti, -,774)6505' 1 = 1 

We now deal with the terms of type 2. By ([20]), if (vri,--- ,114) £ S, ^ULiI^ttiI < 
(^^-g(7ri,7r2,7r3,7r4)+i otherwise (vTi , • • • , 7r4) G S" and, by ([IS]), Enf=i {Ftv^I < Cn-5(''i'''2,7r3,7r4)+2_ 
There are at most Cn'^ quadruples of mappings such that (/(vri, 7r2, vra, 7r4) = q. Hence 

4 

(tti, -,774)6505' 1 = 1 

We turn to the terms of type 3: J2{tti ■■■ TTi)suS'^Fwi^FTri^^FT^i.^Fwi^- Assume first that the 
quadruple (tti,--- ,714) € S'. If h = 1, 'I2 = 3, h = 2, U = 4, then EF^^^F^^^EF^^^F^,^ = 
]^^^^ EFtt; and we obtain the same bound that the terms of type 1. The other cases reduce 
to the case h = 1, I2 = 2, I3 = 3, k = 4 and by ([IB]), EF^^F^^ < C'n-'?(''i'''2)+i_ rj.^^^^ 
are at most Cn'^'^'^ quadruples such that q{TTi,TT2) = q and g(vr3,7r4) = q'. We deduce that 

E{7ri,-,7r4)G5'^^-^^l-^^2lEF^3F^4 < Cn^. 

Assume now that that (vri, • • • , 7r4) G S. We have: 
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If A^(7r;^ , vr^j) = 0, then N{-ki^,'Ki^) > 1 and there are at most Cn'^^'''^'^ quadruples such 
that (/(vr^^, TT/j) = q and (/(vr^gjTr^) = q'. Indeed, since (tti,--- ^tt^) £ S, the cardinal of the 
intersection of the images of (vr/^jVr/j) and (tti^^tti^) is at least 2. The other cases reduce to 
the case, N{tti^,tti^) > 1 and N{iti^,tt^) > 1, for such cases, we notice that there are at most 
(jjiQ+q quadruples such that q{'Ki^,TTi^) = q and q{'Ki^,-Ki^) = q' . In all cases, we conclude that: 
EK,..,.4)eslEi^-i^-2lE^^3^^4 < Cn\ Hence, 

(ni,- ,TT4,)esuS' 

It remains to treat the terms of type 4. Assume that (vri, • • • , 7r4) € S U S", we have: 

KF^,EF^^^EF^^,^F^^^ < C'„-g('rii) + l-g(7ri2,7ri3,7r,^)+e(7r)^ ^^21) 

where e(7r) € {1,2}, e(7r) = 2 if there exists j G {2,3,4} such that A^(7r/^. , vr/^.) = for all 
k G {2,3,4}\{j}, otherwise, e(7r) = 1. 

If e(7r) = 1 then since there are at most Cn'^'^'^ quadruples such that q{TTij^) = q and 
q{iTi2,TTi^,TT^) = q', we deduce that E{7ri,... ,7r4)GSuS' = '^)^Fni^Fni^^F7ri.^Fni^ < Cn^. 

If e(-7r) = 2, then, without loss of generality, we may assume N(TTi^,7ri^) = for k G {3,4}. 
It implies that 5(71^2, vr^g, vr^) = q^ni^) + q{TTi^,TT^). Since (vri, ■ ■ ■ ,114,) £ S U S' , N{tti^,tti^) > 1, 
therefore q{TTi^) + ^(vr/j) > qiT^inT^h) + 1 and by Inequality (f2T]) . 

EF^, EFr, EF^, F^, < Cn"'?('^'i)+^^'?('''2''^i3''^'4)+2 < (-7^-g(7r,,,7r,2)-9(7r,3,7r,J+2_ 

Finally, we notice that there are at most Cn'^^'^ quadruples such that (^(Tr;^, vr/j) = q and 
q{iTi^,iTi^) = q' , it follows that: 

^ EF^,^EF^,^F^,^F^,^ <C7n2. 

(7ri,---,7r4)e5U5' 

Inequality (fTBI) is proved. 

3.2 Proof of Theorem H 

By Equation (|15p . for m > 2, we have (with yi =0): 

m p_i „ m 

i^^{Pm) = f{Or + Y.\- / n/(y^0)-y<^0+i))'^y2---dyp. (22) 

P=2 P- 0es„./(^'')''-S=i 

The leading term in 7 is of order 7*""^. Taking p = m in the above expression gives: 

l^-yiPm) ~ 7"""^ / IT fiVj - Vj+l)dy2 ■ ■ ■ dUm- 

A direct iteration leads to: 

/ n f^y^ - yj+i)dy2 ■ ■ ■ dyrn = /*"^(0), 

where / * g{y) = f{x)g{y - x)dx, /*^(x) = /(x), and for m > 2, = * /. 
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Hence f^^d)m-i lYJli fiVj - yj+i)dy2 ■ ■ ■ dym = f^d f'^iOdi = j t'^i){t)dt and for all m > 2, 

t 



u^iPm) ~ 7™-' J t^mdt = J t^r'i^{-)dt. 

Since / tip{t)dt = J = /(O), this formula is still valid for m = 1. Now, let h(t) = 
Y^m>i ^mt"^ with Ylm>i l^ml*™ finite for all t, then since, \u^{Pm)\ < m7™~-^C™, using Fubini's 
Theorem, the conclusion follows. □ 
Remark. We can easily identify the next term in the asymptotic of i'^{Pm). The second 
leading term in Equation (p2|) is of order 7™"^. As in ([TO]) , it is equal to 

m—l 



I I It ^ li n 

y p vPip{u)du / 
Jr Jm. 



Since if n 7^ w, Y^pLi P^^'w'" ^ = uv{u — v) {{m — l)u^ — mu^ v + v"^), we deduce that: 
3.3 Proof of Proposition [4] 

Let Dn denote the n x n matrix with entry i,j equal to: Idl-'^j — -'^jll < Sn), if In denotes the 
nxn identity matrix, Dn — In is the adjacency matrix of the random geometric graph Q{Xn, dn) 
where there is an edge between i ^ j if \\Xi — Xj\\ < 5n- We have component wise: 

-\\f\\ooDn<B^<\\f\\^Dr,. 

Since the spectral radius p{Bn) of Bn is upper bounded by maxi<j<„ | Yl'^=i{Bn)ij\-, we deduce 
that: 

p{Bn) < ||/||oo(l + A„), 

where A„ is the maximal degree of the graph Q{Xn,5n)- Then, the proposition follows from 
Theorem 6.6 of Penrose [15]. □ 



4 Further properties of the Euchdean Random Matrices 

4.1 Eigenvectors of Euclidean Random Matrices 

As it is pointed by Mezard, Parisi and Zee [13j, if Ui = {^k,n)i = e^'^'"^'-^^ we have: 



JliTvk.x 



In particular, if F{x) = e , then n is an eigenvalue with ^k,n as eigenvector and the 
rank of A is 1. Note also by the Strong Law of Large Numbers that for all i, a.s. 



lim - YF{Xi - x,)e-2-'= (^>-^.) = F{k). 



n 

J 

if An = A/n, by Equation (j23p . for all i, a.s.: 

lim {An^k,n)i = F{k){^k,n)i- 
n—^QO 

This last equation is consistent with Theorem [1) a.s. for n large enough there exists an 
eigenvalue of An close to F{k). It is possible to strengthen this last convergence as follows: 
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Proposition 9 For p > 1, let \\U\\p =( J2i>i\Ui\^ ^'^'^ ll^^lloo = supj>;^ |f7j| . For all 
p G (2, oo], a.s. for all k G Z'^, 

lim - F{k)<^k,n\\p = 0. 

n— +00 

Moreo7;er, lim„^oo EP„$fc,„ - F{k)^k,n\\l = ll/lll " 1^(^)1' = Ei^k 1^(01'- 
Proof. To simplify notation, we write $ = and /jj = F{Xi — Xj). 



n -, n 



\An'^k,n - F{k)<^k,n\\p > £ 



i=i j=i 
1 " 

n 



p eP 
> — 
n 



> en 



i-i/p 



From Equation \ E"=i /ij^i - nF{k)^i\ = \ E"=i ^(^i - Xj)e-^^''''-^^^-^J^ - nF{k)\. 

Hence: 



(^\\An^k,n-F{k)<^k,n\\p>e) < nP ( | ^ - X 



; e 



~2i7rA:.(Xi-Xj) 



< nE 



< 2nexp(- 



I ^ F{Xi - Xj-)e-2^^'=-(^i-^^-) - (n - 1)F(A;)| 
i=2 

> en^-i/p- |F(0)| - |F(A;)| X 

max(0, (en^-i/P - |F(0)| - \F{k)\)f 
lli^lloo(n-l) 



where the last equation is Hoeffding's Inequality. We then apply Borel Cantelli Lemma. 
It remains to prove the statement of the proposition for p = 2. Similarly, we obtain: 



Ep„$fc,„ - F(A:)$fc,n|li = -E| - Xj)e 



'2iTTk.{Xi-Xj) 



nF{kf 



We then write E I Y.jF{Xi-Xj)e-'^^''^-'^^^'^^'i-nF{k)\^ = E[E[| ( F{Xi-Xj)e-'^^-^^-'^^^-^^'^ - 

F{k) ) |2|Xi]] = |F(0) - F{k)\'^ + (n - 1) |F(x)e-2*'^'=-^ - F{k)\^dx. The statement follows. 

□ 



4.2 Correlation of the Eigenvalues 

In this paragraph, we state an elementary lemma on the m-correlation of the eigenvalues of A 
(m < n): 

/ \ 

^- = vrjE y: 11(^^.-^(0))' 

{h,--- ,im}C{l,-n} j=l 

where the sum is over all subsets of {1, • • • ,n} of cardinal m. Note that Mi = and that 
Mm is related to the factorial moment measure pm{dzi, ■ ■ ■ , dzm) (also called the joint intensity 
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measure, refer to Daley and Vere-Jones [6]) of the point process {Ai — F{0), • • • , A,„ — F{0)} as 
follows: 



Mm= / T[ ZjPm{dzi,- ■ ■ ,dZm), 



Heuristically, pmidzi, ■ ■ ■ , dz^) is the infinitesimal probability of having an eigenvalue at F{0) + 
Zi for each i € {!,■ ■ ■ , m}. We define A = A — F{0)I, where I is the nxn identity matrix (note 
that "~" is not the complex conjugate of the matrix A). A{xi, ■ ■ ■ ,Xm) is the m x m matrix 
where the coefficient i,j is equal to F{xi — xj) — dijF{0). 

Lemma 10 

Mm = / detA{xi,--- ,Xm)dxi---dXjn- 

For m = 2 we get: 



M2 = - [ F{xfdx, 



the point process of eigenvalues is thus repulsive. 

Proof. The characteristic polynomial of A is Xa(^) = det(j4 — tl) = YYi=ii\ ~ /(O) ~ t) = 
X]J^=o"m(~*)"~™'' where, = IVjLii^ij - f i^)) ■ However, by Newton formula, we 

also have, = X^jj^ ... det^jj^ ... j^}, where for a set of indices i = {ii, ■ ■ ■ ,im}, A is the 
mxm extracted matrix obtained from A by keeping the raws and columns {ii, • • • , im} (i-e. Ai is 
a principal minor). Taking expectation, we deduce that = (^) /^m det A{xi, • • • , Xm)dxi ■ ■ ■ dxr, 
□ 

In Lemma [TOl we have computed the mean value of the symmetric polynomials: 

m 

{«i,---,«m}C{l, ■■■?!} i=l 

for the vector A = (Ai — F(0), • • • , A^ — -^(0)). Actually, it is possible to compute the mean 
value of the symmetric polynomials: 

m 

^m,k i^i , • • • , ^ ^ Xj . 

{h,--- ,«m}C{l, - n} J=l 

for the vector A. To this end simply consider, Xa*=(0 = det(^'^ — tl) = Y[7=ii(.^i ~ /(O))'^ ~ 0- 
We obtain similarly: 



l/( )Eamfc(A) = / Y\ z'^pmidzi,- ■ ■ ,dzm) = det A'^{xi,- ■■ ,Xn)dxi- ■■ dxn, 



where m = {1, • • • ,m} and for a set of indices i = {zi, • • • ,im}, is the mxm extracted 
matrix obtained from A by keeping the raws and columns {ii, - ■ ■ , im}- 
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